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ABSTRACT 

We describe a new code which can accurately calculate the relativistic effects which 
distort the emission from an accretion disc around a black hole. We compare our results 
for a disk which extends from the innermost stable orbit to 20r^ in both Schwarzschild 
and maximal (a = 0.998) Kerr spacetimes with the two line profile codes which are 
on general release in the XSPEC spectral fitting package. These models generally give 
a very good description of the relativistic smearing of the line for this range of radii. 
However, these models have some limitations. In particular we show that the assumed 
form of the angular emissivity law (limb darkening or brightening) can make significant 
changes to the derived line profile where lightbending is important. This is always the 
case for extreme Kerr spacetimes or high inclination systems, where the observed line 
is produced from a very large range of different emitted angles. In these situations the 
assumed angular emissivity can affect the derived radial emissivity. The line profile is 
not simply determined by the well defined (but numerically difficult) physical effects 
of strong gravity, but is also dependent on the poorly known astrophysics of the disc 
emission. 
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1 INTRODUCTION 

Material in an accretion disk around a black hole is orbiting 
at high velocity, close to the speed of light, in a strong gravi- 
tational potential. Hence its emission is distorted by doppler 
shifts, length contraction, time dilation, gravitational red- 
shift and lightbending. The combined impact of these spe- 
cial and general relati vistic effects was firs t calculated in the 
now seminal paper of ICunninghamI (jlQTST) . where he used a 
transfer function to describe the relativistic effects. The ob- 
served spectrum from an accretion disc around a Kerr black 
hole is the convolution of this with the intrinsic disc contin- 
uum emission. 

While such models have been used to 
try to determine the gravitational potential 
from the o bser ved accr etion disk spectra (e.g 
Laor Net^ Il989l : 
19911: 



■ ■ Ebisa wa. Mitsuda & Hanawa^ 

lEbisawa et all 1199.4 iMakishima et al 2000; 

GierlinskT^Maciolek-Niedzwiecki <§^Ebisawal 1200 lli , these 
attempts suffer from our limited knowledge of the spectral 
shape of t he intrinsic accretion disk emission (see e.g. the 
review It is much easier to determine the 

relativistic effects from a sharp spectral feature, such as the 
iron fluorescence line expected from X-ray illumination of 
an accretion disc (jFabian et al.l[l98^ . An originally narrow 
atomic transition is transformed into broad, skewed profile 
whose shape is given directly by the transfer function. 



Observationally, evidence for a relativistically smeared 
iron line first came from the ASCA observation of the activ e 

I ll99.^ . 



galactic nuclei (AGN) MCG-6-30-15 iTanaka et alJ 
Further observations showed evidence for the line profile be- 
ing so broad as to re quire a maximally spinning black hole 
jlwa sawa enilll99(j) . More recent data from XMM are in- 
terpreted as showing that the line is even wider than ex- 
pected from an extreme Kerr disk, requiring direct extrac- 
tion of the spin energy from the centr al black hole as we ll 
as the immense gravitational potential (jWilms et aLlEoOll) . 

Such results are incredibly exciting, but X- 
ray spectral fitting is not entirely unambiguous. 
There is a co mplex reflected continuum a s well 
as the line jN avakshin, K azanas Kallmlinl l2000l : 
iBallantvne. Ross F abian 200 1^^^ For an ionised disk 
(as inferred for MCG-6-30-15) the current models in general 
use (pexriv in the XSPEC spectral fitting package) ar e 
probably highly incomplete ()Ross. Fabian Yound 119991) . 
Co mplex ionised absorption also affects AGN spectra (e.g. 
.KasDi et al.ll2002^ and the illuminating continuum itself can 
have complex curvature rather than being a simple power 
law. 

However, in MCG-6-30-15 these issues have been exam- 
ined in detai l, and the results on the dramatic line width ai 
pear robust dVaughan Fabianl2003l : lR,evnolds et al.l200 
Thus there is a clear requirement that the extreme relativis- 
tic effects are well modelled. There are two models which 
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are currently widely available to the observational commu- 
nity, withi n the XSPEC spect ral fitting p ackage, dis kline 
(based on iFabian eral] ll 989h and laor (|Laoi1ll99lh . The 
analytic diskline code models the line profile from an ac- 
cretion disc around a Schwarzschild black hole (so of course 
cannot be used to describe the effects in a Kerr geometry). 
Also, it does not in clude the effects of lightbending (although 
iFabian et al.lll989l outline a scheme for incorporating this) 
and hence does not accurately calculate all the relativistic 
effects for r < 20r^ (where rg = GM/c^). By contrast, the 
laor model numerically calculates the line profile includ- 
ing lightbending for an extreme Kerr black hole, but uses a 
rather small set of tabulated transfer functions which limit 
its resolution and accuracy (see Section 3.3). 

While there are other relativistic codes in the literature 
which do not suffer from these limitations, these are not gen- 
erally readily and/or easily available for observers to use. 
There is a clear need for a fast, accurate, high resolution 
code which can be used to fit data from the next genera- 
tion of satellites. In this paper we describe our new code 
for computing the relativistic iron line profile in both the 
Schwarzschild and Kerr metrics. We compare this with the 
diskline and laor models in XSPEC for discs which extend 
down to the last stable orbit in their respective spacetimes, 
and highlight both the strengths and limitations of these 
previous models. 



CALCULATING STRONG GRAVITATIONAL 
EFFECTS 

;. ICunninghamI 1197.4 

and calculate an in- 



We follow the standard approach (e 
IFabian et a'l]ll989l: iFanton et al.lll99 / 

finitesimal amount of flux, dFo observed at energy, Eo due 
to a patch on the disk which subtends a solid angle dE on 
the image of the disc at the observer (see Fig. Q]and|2J. 

dFo {Eo) = lo {Eo) dE = g^Ie {Ee) dE (1) 

where the redshift factor g = Eo/Ee and specific intensity 
in the observers and emitters frame, lo and le are related 
through the relativistic invariant I/u^. For an emission line 
with rest energy Eint, then le {Ee) = £ {re, fie) ^ {Ee - Eint), 
where e (re, /ie) is the emissivity, which can be a function of 
the radius, re and angle, /He at which the photon is emitted 
(as defined in Figure The infinitesimal flux becomes 

4 



dFo {Eo) = g £ {re.lJie) 5 {Eo - gE^nt) dE 



(2) 



The total flux can be obtained by integrating over all the 
entire image of the disk in the observers sky. We can write 
dE = dad(3/rl where a, f3 are the x,y coordinates of the 
image of the disc at the observer with coordinates {ro,0o) 
(see Figure I2J, such that 



Fo {Eo) 



g £ {re.fJie) 5 {Eo - gEint) dad(3 (3) 



The a, p position of the image of the disc section is related 
to the conserved quantities, A,^ which describe the contri- 
butions to the photons angular momentum from the radial, 
polar and azimuthal directions (^de Felice PretL1999') . via: 



A 



P = ±\/q^ COt2 Oo 



(4) 
(5) 




Figure 1. The coordinate system used for the disc. The emission 
is defined in the rest frame of the disc material. The polar and 
azimuthal emission angles are obtained by taking the dot- 

products of the photon four-momentum with the basis vectors of 
this frame, where fie = cosB. This disc frame can be connected 
to the frame which co-rotates with the black hole spacetime via a 
simple boost which depends on the velocity structure of the disc. 




image on 
observers sky 



Figure 2. Diagram showing the link between the observers frame 
of reference and the global coodinate system defined by the black 
hole. Photons that are emitted from the disc at some distance re 
from the hole are seen at coordinates a, (3 on the image of the 
disc at the observer. 



For a thin, Keplerian disc, these constants of motion can be 
written in terms of the redshift factor of the photon, g and 
the radius of emission, re and angle of emission, /ie of the 
photon (as previously defined): 



A 



9 

ref^e 



(6) 
(7) 



Here, Q describes the azimuthal velocity profil e of the emit- 
ting region and is the 'redshift function' (jFanton et al.l 
Il997l : lMartocchiall200(tl . which for a geometrically thin, Ke- 
plerian disc located in the equatorial plane are given by: 



1 



{1-aQ)^ - {rl + a^) 



(8) 



(9) 



Thus the problem reduces to finding the area on the ob- 
servers sky subtended by all parts of the disc which con- 
tribute to a given Eo. 
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Figure 3. Heavy blue lines denote adjacent contours of constant 
redshift, gi < p^+i on the observers sky (the a, /3-plane) that 
define the area of the redshift bin g = -\- ^ {9i-\-i — 9%) with 
width dg = gi-\-i — gi. Red narrow lines show examples of the 
divisions used to create a set of meshed trapezoids that enable 
the area of the redshift bin to be determined. For the purposes of 
clarity, this mesh is far coarser than that used in the calculations. 



The photons (null geodesies) that link the ac- 
cretion disc with the observer can only be found 
by determining the full general relativistic light travel 
paths which link the disc to the observer. These null 
geodesi es ar e given by solutions of t he geodesic equa- 
tions Jc^er 196R: |Misner, Thor^nd Wheeler 1973; 
IChandrasekhai<ll98.^. which can be obtain ed numerically 
(e.g. iKaras. Vokrouhlickv PolnarevI Il992fl. but can also 
be given in terms of analytic funct ions (iRauch BlandfoTdl 
ll994l : lAgo]|[l997l : ICadez. Fanton Rr, Calvanilll998^ . which en- 
able them to be solved quickly and with arbitrary accuracy. 

We develop a technique based solely on the image 
of the accretion disc in the a,p plane, which defines the 
flux received by the observer, si milar to that employed by 
ICadez. Calvani FantonI (|2003f) . This allows us to gener- 
ate high resolution, accurate line profiles numerically while 
avoiding the issues surr ounding the pa rtial derivatives of 
the geodesic equations (IViergutj[l99^ . We use the ana- 
lytic solutions of the geo desic equations as tabulated by 
iRauch BlandfordI (|l994r) to find the complete set of light 
travel paths that link the accretion disk and the observer at 
(vo^Oo). We sort these by redshift factor, and use adaptive 
griding to find the boundarys on the (a, ^5) -plane for all lines 
of constant g. 

Two adjacent boundaries, gi and gi+i, therefore define 
the area of the redshift bin g = gi-\- ^ (gi+i — gi) with width 
dg = gi+i — gi when projected onto the (a, (3)- plane (as is 
shown in Figure E| • We can simply determine the area of 
this region by dividing it up into a set of tessellating trape- 
zoids, as shown in Figure the area of of each of which can 
be determined by a simple geometric formula. The final area 
of the redshift bin is determined by summing together the 
contributions from all such trapezoids internal to {gi^gi+i). 



Each individual trapezoid is small, so that there is no signif- 
icant change in Ve or /ie (though this is not necessarily true 
across the total area dadf3). The emissivity law can be con- 
volved into the calculation using the emission coordinates 
at the centre of each trapezoid to weight its area before per- 
forming the summation over all trapeziods. This approach 
allows us to calculate line profiles at high spectral resolution 
on timescales of a few minutes on a 2GHz desktop PC. 

We have extensively tested the routines that calcu- 
late t he null geo desic paths against those supplied by Eric 
Agol (|Ago]|ll99^ and have found them to be indistinguish- 
able. We have also compared the line profiles generated by 
our code to those presented previously in the literature, 
in particular those generated from the code described by 
iFanton et alJ (|l997l 'l and have again found them to be indis- 
tinguishable. 



3 HIGHLY RELATIVISTIC LINE PROFILES 
3.1 Introduction 

We have taken a disc from rms to rmax — 20r^ (beyond 
which strong gravitational effects become of diminishing im- 
portance) for both the Schwarzschild (a = 0, rms — 6r^) 
and maximal Kerr (a = 0.998, rms = 1.235r^) cases for 
Oo = 5°, 30°, 60° and 85°. 

We first consider the extent of lightbending effects 
in a Schwarzschild spacetime. Figure 3] shows the three- 
dimensional surface in (/Xe, Te, g) for the complete set of light 
travel paths connecting the accretion disc to the observer. 
There is a considerable range of /ie contributing to the ob- 
served emission at all inclinations. For low inclinations the 
effect is fairly uniform, with each radius contributing a simi- 
lar range in /ie, but with a systematic shift to larger emission 
angles (smaller fie) with smaller radii. By contrast, at higher 
inclinations the lightbending is strong enough to gravitation- 
ally lens the far side of the disc. This leads to a much larger 
range of /He which contribute to the disc image at small radii. 
In all cases, lightbending means there is a range of fie which 
contribute to the observed disc emission, so that in general, 
the fine profile will depend on the angular distribution of 
the emitted flux. 

Fig.|5]shows the corresponding surfaces for the extreme 
Kerr case. The disc now extends down to 1.235r^, far closer 
to the corresponding event horizon than in the Schwarzschild 
case. This introduces a greater complexity to the geodesic 
surfaces. The range of emission angles is from zero to unity 
in all cases, including the nearly face on disc at 5°, which 
has important consequences for the calculation of the line 
proflle. 

To construct the relativistic line proflle, we map these 
surfaces on to the (a,/3) plane as discussed in the previous 
section, forming images of the accretion disc, as have 
been p re viously ca lc ulated bv e.^. [ Cunnin gham &: BardeenI 
(Il97.t): iLuminetJ (ll979^: iHollvwood Melial (ll997D 
'Fant on et alJ (Il99/t l: iFalcke. Melia. Agoll (l200(tl . In Fig- 
ure El we present images of the Schwarzchild disc for the 30° 
(top row) and Oo = 85° (bottom row) cases (others provide 
no new information qualitatively). Images on the left-hand 
side of the flgure are coloured by values of the redshift 
factor, g, as deflned by the scale at the top of each image. 
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r, 



Figure 4. Surfaces in the {^e,re,g) parameter space describing geodesies that reach an observer at a given inchnation for a standard 
accretion disk around a Schwarzschild black hole. Notice that, for every (re, g) pair, there are two geodesies that reach any given observer, 
corresponding respectively to geodesies that are emitted from the side of the disc closest to the observer (lower surface) and geodesies that 
are emitted from the opposite side of the disc to the observer (upper surface). These are the two geodesies referred to by Cunningham 
(1975). 




Figure 5. As in Figure^ for a maximal (a = 0.998) Kerr black hole. These surfaces are exhibit a far greater complexity than those in 
the Schwarzschild case. The range of accessible redshift is increased for a given disc-observer system, whilst the range of emission angles 
is from zero to unity for all inclinations. 
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Figure 6. Redshift images (left hand panels) and flux image (right hand panels)of the accretion disc on the (a, (3) plane for a Schwarzschild 
black hole, for the 60 = 30° (top row) and 60 = 85° (bottom row) cases. Redshift images are coloured by the associated values of g 
as measured by the distant observer. Flux images colored by the g^r'^dE. = g^dad(3 component of the relativistic line proflle. Note the 
appearance of strong light bending effects in the Oo = 85° case, as previously reported by Matt, Perola and Stella (1993), Zakharov and 
Repin (2003). 



By contrast, images on the right-hand side of the figure 
has each redshift bin coloured by its area on the observers 
sky, i.e. g^dadf3, again with the scale defined at the top of 
each image. Strong gravitational lensing effects can now 
clearly be seen in the high inclination images. Photons 
from the far side of the disc pass clos e to the black hole, so 
the disc image is strongly di^ orted (lMa,tt. Perola S,^, Stellal 
ll993l : IZakharov ReDinll2003il . Since the area of the disc is 
magnified, its contribution to the observed flux should be 



large. However, we stress that the low inclination images 
are also affected by lightbending (see Figure EJi and^t), 
though they are not magnified by gravitational lensing. 



The form of the line profile is now determined from 
the flux image (representing the effects of strong gravity), 
together with the assumed form for the emissivity (deter- 
mined by the energy release and radiative transfer pro- 
cesses), which is generally taken as (ignoring azimuthal de- 
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Redshift, g Redshift, g 

Figure 7. Comparison of the relativistic hne profile computed by our model (red solid line) with that computed by the XSPEC diskline 
model (blue dashed line) for e (re) oc and / (/ie) = 1- At inclinations of < 30°, the profiles match to within ~ 10%, but the increaing 
importance of lightbending (which is not included in the diskline code) gives a 40% discrepancy in the profile shapes for inclinations 
> 60°. In this and all subsequent figures the line profiles are normalised such that they contain one photon, and all our results are 
unsmoothed. 



pendence) : 

£(re,/ie) = e(re)/(/ie) (10) 

While the flux image is a difficult numerical problem, it 
depends on well known physics. By contrast, the emissiv- 
ity laws considered have rather simple forms, but are de- 
termined by the poorly known astrophysics of the disc. Of 
course, there are many other outstanding theoretical issues 
that can produce a substantial impact on the line profile, 
including (but not Umited to) returning radiation or light- 
bending th at can enhance the emissivity of the inner part 
of the disc (ICunninghamlll975l: iLaor. Netzer PiranlHoOol: 
iMartocchia. Kara'^^Mattll200(tl . emission from the plung- 



ing region (|Revnolds Bege lman||l 997r) and azimuthal de- 
pen dence of the emissivity tead ez. Calvani FantonlEoo3 : 
Kar as. Martocc hia & Subr 2001). However, these are out- 
side the scope of the current work. 



3.2 Comparison with the Diskline Model 

The diskline code assumes a Schwarzchild metric (a = 0) 
and additionally that light travels in straight lines (so the 
angular emissivity term is irrelevant). In its XSPEC imple- 
mentation it allows both arbitrary power law e (re) oc 
and point source illumination. However, its analytic struc- 
ture means that any radial emissivity law is easy to incorpo- 
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Redsliift, g 

Figure 8. Comparison of the relativistic line profile computed by 
model (blue dashed line) for e (re) oc and / (/ie) oc (1 + 2.06/ie)- 



rate. We choose to use q — —3, as this is approximately the 
form of the gravitational energy re lease per unit disc area 
(see e.g.lZycki, Done PmithjlHi). 

Figure EJ shows our line profiles assuming /(/ie) = 1 
(no angular dependance of the emissivity) compared with 
those from the diskline code. We see that our new model 
matches very closely to the XSPEC diskline model for a 
nearly face on disk. Whilst the key difference between our 
model and diskline is the inclusion of light-bending effects, 
the impact of this is small at low inclinations if there is no 
angular dependance to the emissivity (but see Section 3.4). 

By contrast, at high inclinations, lightbending not only 
means that the line is formed from many different fie, but 
gravitational lensing enhances the flux from the far side of 
the disc. This lensing effect gives clear differences between 
our model and diskline. The lensing magnifies the image of 




0.0 0.5 LO L5 

Redshift, g 



our model (red solid line) with that computed by the XSPEC laor 
The profiles produced by the two models match to within 5 — 10%. 



the far side of the disc, which has velocity mostly tangential 
to the line of sight, so is not strongl y doppler shifted. This 
boosts the line profile at ^ 1 (see iMatt. Perola Stellal 
1993). Since the line profiles are all normalised to a single 
photon, then this also makes the blue peak smaller. 

In summary, the diskline model as incorporated into 
XSPEC produces line profiles which are accurate to ^ 10% 
for inclinations of less than 30° . Obviously, if the inner disc 
edge rmin > Tms then the lightbending effects become cor- 
respondingly smaller and the match between the two codes 
becomes even closer. At higher inclinations the differences 
between diskline and our code become larger due to the 
effects of gravitational lensing, which leads to an effective 
redistribution in flux between the blue peak and the center 
of the line compared to that predicted from straight light 
travel paths. 
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Figure 9. Comparison of the relativistic line profiles generated 
by our model with (a) e(re) oc r^^, /(^e) = 1 (red lines), (b) 
e (re) oc r~^, / (/ie) oc (blue lines), (c) e (re) oc r~^, / (/ie) oc 
(1 + 2.06/ie) ( lines) for a maximal Kerr black hole with 

the disc extending from 1.235 — 20rg. The relative height of the 
blue wing changes by ~ 35% for different angular emissivity laws, 
anti-correlated with the slope of the red wing. 



3.3 Comparison with the Laor Model 

By contrast, the laor code is written for extreme Kerr, and 
includes a standard limb darkening law /(/ie) oc (1 + 2.06/ie). 
The code is based on a series of photon trajectory calcula- 
tions, where the disk is split up into a set of rings of width 
dre at re . Each part of the ring radiates with total emissiv- 
ity (radial plus angular) given simply by the limb darkening 
law (i.e. no radial dependance, ^ = 0) and the Une profile 
from that ring is built up from many light travel paths which 
connect the disc to the observer. This produces a series of 
transfer functions T{re,Eo — gEint) at each radius, analo- 
gous to Figure ISJi-d but including the limb darkening law. 
These tabulated transfer functions are read by the laor code 
in XSPEC and used to build a total line profile for any given 
radial emissivity Fo{Eo) = J e(re)T(re,Eo - gEint)redredg. 

We compare this with our code, using a ^ = — 3 emis- 
sivity for both as in the diskline comparisons above. We 
include the same limb darkening law as used by laor and 
the results (Figure |H| show that the overall match between 
our code and laor is good to 5 — 10%. 



4 THE ROLE OF THE ANGULAR 

EMISSIVITY AND BLACK HOLE SPIN 

The effect of appling a radial emissivity is straightforward. 
The transfer function describing all the relativistic effects 
from a given radial ring of the disc is unaffected, so the 
effect is simply to change the weighting of the line profile 
from each radial ring of the disc. 

By contrast, the effect of the angular distribution is far 
more subtle. A given radial ring on the disc can contribute to 
the line profile from a range of emission angles. The relative 
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Figure 10. As in Fig.|n]witli the disc now extending from 1.235 — 
400rg. There is still a ~ 25% change in the blue wing height 
and significant change in red wing slope for the different angular 
emissivities, despite the inclusion of the outer disc regions. 




Redshift, g 

Figure 11. As in Fig.lolwith the disc now extending from 1.235 — 
6rg. The additional magenta line is for a limb darkened angular 
emissivity with more centrally concentrated radial emissivity, oc 
r^^'^ . This is very similar to the blue line profile derived from 
a very different radial emissivity, oc r^^ , with a limb brightened 
angular emissivity. 



weighting of these is determined by the angular emissivity, 
so it forms part of the calculation of the transfer function 
itself. 

Different angular emissivity laws can have striking ef- 
fects on the form of the relativistic line profile, which we 
illustrate in Fig. |^ for a maximal Kerr geometry (a=0.998) 
with the disc extending (as previously) from 1.235 — 20r^ 
and 6o = 30°. The line profiles here all implement the stan- 
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Figure 12. As in Fig.Elbut with the disc extending from 6— 400r^ 
in an extreme Kerr (sohd hne) and Schwarzchild (dashed hne) 
spacetime. The differences between the hne profiles produced for 
the same sized disc in different assumed spacetimes is of order 
~ 5% for a given angular emissivity. The effect of changing the 
angular emissivity is also similarly small (~ 5 — 10%). This con- 
trasts with the much larger effects seen in the extreme Kerr metric 
for a disc extending down to 1.235rp, where lightbending is much 
more important (Fig. EJ. 




0.7 0.8 0.9 

Redshift, g 

Figure 13. As in Fig. 1121 but with the disc extending from i 
20r^. The differences now are of order 15-20%. 



dard radial emissivity law of r~^. However, we now compare 
a range of angular emissivity laws, these being (from top to 
bottom at the blue peak in Fig. |3 the standard limb dark- 
ening law (as discussed in §3.3), followed by the constant 
angular emissivity case (as used in §3.2). An ionized disc 
could also be limb brightened^ with the probable limiting 



case of /(/i) cx: l//i as expected from optically thin material, 
shown as the bottom line in Fig. 1^1 There is a 35% differ- 
ence in the height of the blue peak depending of the form of 
the angular emissivity used. 

However, such a limited range of radii is probably not 
very realistic. The disc should extend out to much greater 
distances from the black hole, where the relativist ic effects 
(including lightbending) are less extreme. However, realistic 
emissivities strongly weight the contribution from the in- 
nermost regions, so the effective dilution of the relativistic 
effects by including the outer disc is not overwhelming. Fig. 
QUI shows the line profiles generated using the same angular 
emissivity laws for a disc extending from 1.235 — 400rp , again 
with Oo — 30°. There are still significant differences in the 
line profiles, with a ^ 25% difference in the height of the blue 
peak while the red wing slope changes from Fo{Eo) oc E^'^ 
(Umb darkened) to oc E^'^ (Umb brightened). 

Despite the expectation of an extended disc, som e re- 
cent observational studies (e.g. iR.evnolds et al.l 12004 ) have 
tentatively suggested that the disc is very small, from ~ 
1.235 — 6rg. This enhances the importance of lightbending. 
Fig. 111! shows the line profiles for a disc extending from 
1.235-6r^ using the different angular emissivity laws of Fig. 
El The blue peak height differences are ^ 40%, and the red 
wing slopes are different. For comparison we also show a 
Umb darkened profile obtained from a very different radial 
emissivity of r~^'^ . This is very similar to the extreme limb 
brightened profile obtained from the r~^ radial weighting. 
We caution that uncertainties in the angular distribution of 
the line emissivity can change the expected line profile due 
to lightbending effects even at low/moderate inclinations, 
and that this can affect the derived radial emissivity. 

Currently, the only available models in XSPEC have ei- 
ther zero or maximal spin. A zeroth order approximation to 
spacetimes with different spins is to use the maximal Kerr 
results but with a disc with inner radius given by the min i- 
mum stable orbit for the required value of a fe.g. lLaoJlQQlh . 
We test this for the most extreme case of a = modelled 
by a maximal Kerr spacetime with rmin = 6r^. Fig.lT^com- 
pares this with a true Schwarzschild calculation for a disc 
extending from 6 — 400r^ with Oo = 30° for a range of an- 
gular emissivities. The differences between the spacetimes 
(for a given angular emissivity) are at most 5%. This is 
roughly on the same order as the effect of changing the an- 
gular emissivity, which is much reduced here compared to 
Fig. 1^ due to the larger rmin- Assumptions about both spin 
and angular emissivity become somewhat more important 
for smaller outer disc radii. Fig. shows this for a disc 
between 6 — 20r^ (directly comparable to Fig. [7J. 



5 CONCLUSIONS 

Recent observational studies have provided evidence for 
highly broadened fluorescent iron K a lines. While there 
are a variety of line profiles seen (e.g. iLubinski Zdziarskil 
I2OO J) , there are some objects where the line implies that 
there is material down to the last stable orbit in a maxi- 
mally spinning Ke rr spacetime (most notably MCG-6-30-15: 
IWilms et aTll200lh . However, the strong gravity codes gen- 
erally used to model these effects are now over a decade old. 
Increased computer power means that it is now possible to 
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improve on these models. We describe our new code to cal- 
culate these effects, which uses uses fully adaptive gridding 
to map the image of the disc at the observer using the ana- 
lytic solutions of the light travel paths. This is a very general 
approach, so the code can easily be modified to incorporate 
different emission geometries. 

We compare the results of our new code with those 
from diskline and laor (publically available in the XSPEC 
spectral fitting package) for Schwarzchild and extreme Kerr 
spacetimes. These previous models are accurate to ^ 10% 
with realistic (oc r~^) radial emissivities. However, they 
make specific assumptions regarding the angular dependence 
of the emitted flux, which may or may not be valid. Light- 
bending is always important for a disc which extends down 
below 20r^, in that the image of the disc at the observer 
always consists of a range of different emission angles. This 
can produce significant changes to the derived line profile, 
especially in extreme Kerr spacetimes. Whilst calculating 
strong gravitational effects is a difficult numerical problem, 
the underlying physics is well known. By contrast, the angu- 
lar emissivity is an astrophysical problem, and is not at all 
well known as it depends on the ionization state of the disc 
as a function both of height and radius. Before we can use 
the line profiles to provide a sensitive test General Relativ- 
ity and probe the underlying physics, we will need to have a 
much better understanding of the astrophysics of accretion. 

This code will be publically released for inclusion as a 
convolution model in the XSPEC spectral fitting package. This 
will include arbitrary spin and inner and outer disc radii as 
well as allowing both angular and radial emissivities to be 
specified. 

After this paper was submitted we learnt of the indepen- 
dent work by (jPovciak. Karas YaaoobI 120041) which also 
develops a new strong gravity code. Their results match very 
closely with those presented here. 
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